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<D . —o 
c~| , We study semi-inclusive charmless decays B — » ttX in detail, such as B 



X 



ir^X, B° -> vr ±(0) X, B ± -> vr ±(0) X, where X does not contain a charm 
(anti) quark. We find that the process W -> tt-X {B° -» ir+X) can be 
particularly useful for determination of the CKM matrix element \V u b\- We 
calculate and present the branching ratio (BR) of B~° 7r~X as a function of 
\V u b\, with an estimate of possible uncertainties. It is expected that the BR is 
an order of 10 -4 . Our estimation indicates that one can phenomenologically 
determine \V u b\ with reasonable accuracy by measuring the BR ofW^ tt-X 



(B° -> vr+X). 
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I. INTRODUCTION 



The source of CP violation in the Standard Model (SM) with three generations is a 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix (TJ. A precise measurement of the 
CKM matrix elements is one of the key issues in the study of B— mesons and B— factory 
experiments. In particular, the accurate determination of V u b is one of the most challenging 
problems in B physics. Its non-vanishing value is a necessary condition for CP violation 
to occur in the SM and its accurate value can put strong constraints even on the unitarity 
triangle: for instance, on the magnitude of the CP violating phase (3{= <f>\). 

Theoretical and experimental studies for probing V u b have been mostly focused on the 
semileptonic B— meson decays. The present best experimental data for V u b come from 
measurements of the exclusive decay B — > plv and the inclusive decay B — > X U W, but these 
measurements suffer from large uncertainties due to model-dependence and other theoretical 
errors. For example, the CLEO result obtained using the exclusive semileptonic decay 
B - pW [| : 

\V ub \ = (3.25 ± 0.14(stat.)+°;^(syst.) ± 0.55(model)) x 1(T 3 . (1) 

The OPAL data obtained using the inclusive decay B —>■ X u lv || : 

\y uh \ = (4.00 ± 0.65(stat.)±°;^(syst.) ± 0.19(HQE)) x 10~ 3 . (2) 

The method using the exclusive semileptonic decays involves hadronic form factors, such as 
pB~>-K Qr j[B~>p^ w h oge va l U es are heavily model-dependent and cause large uncertainties. 
The difficulty of using the inclusive charmless decay B — > X U W is in discriminating this 
process from the dominant background B — > XJu decay 0], whose branching ratio is more 
than 50 times larger than that of B — > X U W. 

Although traditional difficulties with the understanding of non-leptonic B decays have 
prevented their use in determination of the CKM matrix elements, the possibility of mea- 
suring \V u b\ via non-leptonic decays of B— mesons to exclusive or inclusive final states has 
been also theoretically explored J5|-[£|. 
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In this work we study semi-inclusive charmless decays B — > txX and investigate the 
possibility of extracting |V^| from these processes. Compared to the exclusive decays, these 
semi-inclusive decays are generally expected to have less hadronic uncertainties and larger 
branching ratios. There are several possible processes in B — > ixX type decays, such as B — ► 
7]-±(°)X, -B — > tt ±( -°^X, B ± — > ir^^X, where X does not contain a charm (anti)quark. The 
class of nonleptonic two-body decay modes and their advantages within general arguments 
were previously discussed in Ref. ||, and semi-inclusive two-body decays within the QCD 



factorization were also studied in Ref. 10 



In Sec. II we first classify all those B — > nX processes, and we identify a certain mode, 
B° — > n~X, whose analysis is theoretically clean and which can be used for determining 
\V u b\. Then, in Sec. Ill we study the mode B° — > tt~X in detail and propose a method 
to extract \V u b\. That is, we calculate the branching ratio (BR) of this mode using the 
full effective Hamiltonian in the framework of the generalized factorization, and present the 
result as a function of \V U b\ with an estimation of possible uncertainties. We also consider 
the B° — B° mixing effect through B° — ► B° — > n~X. The conclusions are in Sec. IV. 



II. CLASSIFICATION OF SEMI-INCLUSIVE CHARMLESS B -» irX DECAYS 

Among the semi-inclusive charmless B — > ttX decays, let us first consider the mode 
B° —>■ n~X. Contributions for the decay amplitude of this mode arise from the color- 
favored tree (b — ► uud) diagram and the b — > d penguin diagram (see Fig. 1), and the tree 
diagram contribution dominates. The charged pion tt~ in the final state can be produced 
via W— boson emission at tree level and is expected to be energetic (E n - ~ m B /2). The 
decay amplitude can be approximated as 

A(B° -> 7r~X) ~ A(b -> 7r~u) ■ h(ud -> X{ud)) w A(6 -> 7r _ «) , (3) 

where h denotes a hadronization function describing the combination of the ud pair to make 
the final state X. To obtain the decay rate, X(ud) should be summed over all the possible 
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(a) (b) 

FIG. 1. Feynman diagrams of 2?° — > ir~X decay: (a) the color-favored tree diagram, and (b) 
the b — * d penguin diagram. 

states, such as 7r + 7r°, 7r + 7r + 7r~ etc, so this process is effectively a two-body decay process 
of b — > n~u in the parton model approximation^ Thus, in this specific mode, no hadronic 
form factors (except the pion decay constant /„-) are involved, and as a result the model- 
dependence does not appear to be severe. We note that the energetic charged pionQ 7r~ in 
the final state can be a characteristic signal for this mode. We will show in Fig. 3 the decay 
distribution, for B° — ► n~X as a function of the charged pion energy E n . (For a detailed 
explanation, see Sec. III.) Like a two-body decay, a peak appears around E n - = ms/2. 

Now let us consider the mode B~ — > it~X. As shown in Fig. 2, various contributions 
are responsible for this process: the color-favored tree diagram, the color-suppressed tree 
diagram, the b — > d and b — > s penguin diagrams. The color-favored tree contribution (Fig. 
2(a)) and b —>■ d penguin contributions (Fig. 2(c)) are similar to those in B° — ► n~X, which 



1 We notice that the dominant tree contribution (Fig. 1(a)) is diagramatically similar to the 
inclusive semileptonic decay, b — >• [l~~v\u, and can be approximated to the free quark decay of 
b — > ir~u within the HQET. 

2 The net electric charge of X should be positive so that such energetic tt~ cannot be produced 
from the inclusive X = 7r + 7r , ir + ir + ir~ , etc. 
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(c) (d) 



FIG. 2. Feynman diagrams of B~ — > ir~ X decay: (a) the color-favored tree diagram, (b) the 
color-suppressed tree diagram, (c) the b — > d penguin diagram, and (d) the b — > d and b — > s 
penguin diagram. 

are effectively two-body type (6 — > 7r~-u) processes. However, the color-suppressed tree (Fig. 
2(b)) and other penguins (Fig. 2(d)) differ from those in B° — > tt~X. In fact, these two 
diagrams correspond to effectively a three-body decay process of B~ — > n~uu (or n~dd) in 
the parton model approximation, with the decay amplitude 

A(B~ -> vr~X) ~ A(5~ -> tt^mm) • -> X{uu)) « A(5~ -> 7r~tm) . (4) 

In Figs. 2(b) and 2(d), the charged pion tt~ in the final state contains the spectator antiquark 
u. So the analysis involves the hadronic form factor for the B — > 7r transition which is 
highly model-dependent. Furthermore, the b — > s penguin contribution (Fig. 2(d)) is not 
suppressed compared to the tree contributions, but dominant in this mode. The decay 
distribution, for £>~ — > 7r~X will be shown in Fig. 4, as a function of the charged pion 
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energy E n . The three-body type contribution from the b — > s penguin is the dominant one. 
Therefore, compared to the case of B° — > n~X, the analysis of this mode is much more 
complicated and involves larger uncertainties. 

Other modes of the type B — > nX can be similarly classified. For instance, in the mode 
B° — > 7r~X, the color-favored tree (b — > uud and 6 — > -u-us) diagrams and the b — > and 
6 — > s penguin diagrams are responsible for the decay process. In this case, the charged 
pion 7r~ contains the spectator quark d so that the process is effectively a three-body decay 
B° — > 7t _ mc? (s) and the hadronic form factor for the B — > 7r transition is involved. Other 
processes are basically a combination of the two-body decay process {b — > 7rg) and the 
three-body decay process — > n~qq'). 

III. ANALYSIS OF if -► vr X DECAY 

In the previous section, we have seen that the process U° — > 7r~X is particularly inter- 
esting, because it is effectively the two-body decay process b — > 7r~w in the parton model 
approximation, and no uncertainty from hadronic form factors is involved. Thus, its theo- 
retical analysis is expected to be quite clean. 

The relevant AB = 1 effective Hamiltonian for hadronic B decays can be written as 



m 



eff 



= 5| [v ub V: q ( Cl O q lu + c 2 0\ u ) + V cb V: g ( Cl Ot + c 2 Ol) 

10 

- E + v cb v^ + v tb v*4) o\ 

i=3 

+ H.C. , (5) 



where OVs are defined as 



0\ f = q^Lff^Lb, 0\ f = WpLfpfrfLba , 
O q 3{5) = qi.Lb E q'YHRW, Ofa = q a ^Lb p E ^L{R)q' a , 

q> q> 

O q m = lqi»LbJ2e q >qYR(L)q', O 8 9 {10) = ^q a ^Lb p J2 e ^ R ( L Wa , (6) 
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where L(R) = (1 =F 7s), / can be u or c quark, q can be d or s quark, and q' is summed over 
u, d, s, and c quarks, a and (3 are the color indices, q's are the Wilson coefficients (WC's), 



and we use the effective WC's for the process b —> dqq' from Ref. [IT]]. The regularization 
scale is taken to be fi = nif,. The operators 0\, 2 are the tree level and QCD corrected 
operators, 3 _ 6 are the gluon induced strong penguin operators, and finally O 7 _i are the 
electroweak penguin operators due to 7 and Z exchange, and the box diagrams at loop level. 

Now we calculate the decay amplitude for the semi- inclusive decay B° — > n~X, where 
X can contain an up quark and a down antiquark. In the generalized factorization approx- 
imation, the decay amplitude is given by 

M = (7i-X\H eff \B°) , (7) 
= i^U{X\u [r(l + 7 5 ) + 1(1 ~ 7 5 )] b\B°) , (8) 

where we have defined the followings: 

r = rribWi , 



2 

mi 



I = -m u w 1 w 2 , 

m d + m u 

Wl = V ub V: d (^- + c 2 ) + ^ + A 4 + ^ + A 10 , 
W2 = -2(^ + A, + ^ + A,) . (9) 



.N c N L 
Here A^ r denotes the effective number of color and 



a = - £ ^ 6 y; d cf. (10) 



We have used the relations 



(7r-|d 7 V«|0> = , (11) 

(TT-i^lO) = -t /<r ™* , (12) 

where and are the decay constant and the momentum of pion, respectively, and m n (m;) 
is the mass of pion (i quark). In Eq. ( |12"D the free quark equation of motion has been used. 



Then, 
\M\ 2 



x 



f^(X\u r(l + 7 5 ) + /(l- 7 5 ) b\W) (2^5\p B -p^-p x ) 



G 



F fl^ B °\ J \X){X\^\B u ){2nY5\ PB -p^ Px ) 



i-prO 



(13) 



x 



where 



r(l + 7 5 ) +/(l-7 5 ) 



b . 



(14) 



In the parton model approximation we take the leading order term in the product of the 
above matrix elements, which corresponds to interpretation of the above process as b(p b ) — > 
^"(Ptt) + u(p u ) Then, \A4\ 2 can be expressed as 



where 



and 



\M\ 2 = 2G 2 F f 2 [(|r| 2 + \l\ 2 )(p u ■ p b ) + 2Re{rl*)m u m b 
= \V ub \ 2 M 2 +\V ub \M 1 + M , 



M 2 = 2G 2 F fl\V* d \ 2 {\c ut \ 2 X + A\c ut \ 2 y + 2Re[(c ut )(c ut y]Z} , 

M 1 = 2G 2 F f 2 \V ud V cb V: d \{2Re[e t ^c ct )(c u y]X + 8 J Re[e^(c ct )(5"')1^ 

+ 2Re[e i \c ct ){c ut Y + e~ l \c ut )(c ct y]Z} , 
M = 2G 2 F f 2 1 V cb V* d \ 2 {\c ct \ 2 X + 4\c ct \ 2 y + 2Re[(c ct )(d ct y]Z} , 



X = (m b + m u )(p b ■ p u ) - 2m b m u 



y 



z 



, tit 



ml 



;(Pb-Pv 



(m d + m u f 

2m u m\ 2 
; {Pb-Pu-m b ) 



c 



N r 



( C 9 C 9) ( C 10 C lo) > 



1 ( 1 

— ^y~( c 3 — C 3) — ( C 4 — c 4) ~~ jv"^ 9 ~~ _ ^ Cln _ Cl 



'10 ) 



= ]f{4 - 4) + (4 - 4) + ^(4 - 4) + (eg - 4) (? = «, 4 . 



(15) 



(16) 



(17) 



Here we have used the usual definition of the phase angle, 

7 (= 3 ) = Arg[-(V ud V: b )/(V cd V; h )] 
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FIG. 3. ^- (in units of 10" 16 ) versus E n for B — > 7r X decay. The solid, dotted, and dashed 
lines correspond to -/V c = 2, 3, 5, respectively. 



We first calculate the decay distribution in the b quark rest frame and boost it to the 
B— meson rest frame. In the b quark rest frame, the decay distribution is given by 

<n 1 M'^—Mjnu /•:, /•:„). (18) 



dE w 



16n rn b E u 



In the B— meson rest frame, the b quark is in motion and has the energy Ej, satisfying 



the relation: E b = m B — E u , where Ei = ^Jp 2 + mf. Ei and rrii denote the energy and 
the mass of i quark inside the B— meson, respectively. The 3-momentum p is defined by 
P — \Pb\ = \Pu\- Thus, the b quark mass is now a function of p given by 



(p) = m| + ml - 2m B \Jp 2 + ml . (19) 



8.0 



Nc-2 

^3body_ds Nc=3 

Nc=5 




0.0 1.0 2.0 



EjGeV) 

FIG. 4. ^- (in units of HT 16 ) versus E n for B — > 7r X decay. The solid, dotted, and dashed 
lines correspond to N c = 2, 3, 5, respectively. Here, '2body_u' stands for the two-body type 
(b — > tt^u) processes (see Fig. 2(a,c)), and '3body_ds' stands for the three-body type processes 
from the color-suppressed tree and the b — > d and b — ► s penguins (Fig. 2(b,d)), while '3body_d' 
stands for the three-body type process from the b — > d penguin only (Fig. 2(d)). 



The decay distribution in the B rest frame can be calculated by 



dT 



Pmax dT^ 

dp p <p(p) 



(20) 



dE n Jo dEn 

where p max — { m %—' m u)/('2 m B)- We have used the ACCMM model [0] using the B— meson 
wave function: 



<f>(p) 



7Cp 3 F 



-_ e ~P 2 /P 2 F 



(21) 



with the normalization J °° p > cj)(p)dp = 1. The Fermi momentum p F = 0.3GeV has been 
used. In Figs. 3 and 4 we show the decay distributions for B° — > n~X and B~ — > ix~X as a 
function of the charged pion energy E n . For B~° — > tt~X decay, the distribution has a peak 
at E n ~ i71b/2. This is a characteristic of a two-body decay. Because the b quark inside 
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TABLE I. The branching ratio (BR) of B -» ?T X for fixed 7 = 65°. Here B denotes the BR 
calculated for \V u b\ = 0.004. 



N c 


B x 10 4 


B 2 x 10 4 


Bi x 10 4 


B x 10 4 


2 


1.05 


0.89 


0.14 


0.016 


3 


1.17 


0.98 


0.17 


0.019 


4 


1.24 


1.03 


0.19 


0.020 


5 


1.27 


1.06 


0.19 


0.021 



the B is in motion, the distribution has some width shown in Fig. 3. For B~ — > n~X, the 
decay distribution is a mixture of three-body type decay distributions and a two-body type 
decay distribution. Figure 4 shows that the dominant contribution arises from a three-body 
type decay (the b — > s penguin process), as explained in Sec. II. 
The BR of 5° — > ix~X can be expressed as a polynomial of \V u b\: 



B(B° -> tj-X) 





2 


B 2 + 




0.004 




0.004 



Bt+Bo, (22) 



where for convenience we have scaled \V u b\ by the factor 0.004 (the central value of the 
OPAL data). Tables I and II show the BR of 5° — > ix~X for fixed j(= ^3) = 65° and 
N c = 3, respectively, with a fixed input value of \V u b\ = 0.004. We note that the term B 2 is 
the dominant contribution (~ 10~ 4 ) to B, while the contribution from the term Bq is very 
small (~ 10~ 6 ). This is due to the fact (see Fig. 1) that B 2 corresponds to mostly the tree 
contribution (b — > uud), while Bq corresponds to the pure b — » d penguin contribution which 
is very small compared to the tree contribution, and B\ corresponds to the interference 
between them. We see that the BR of — > ix~X is about 10~ 4 for different values of 
7(= (f> 3 ) and N c . 

In Fig. 5, we present the BR of — > n~X as a function of N c for three different values 
of 7(= (ps) = 65°, 85°, 110°. As one can see from Eqs. ( |T5| , |1~6| , [22] ) and from Table 2 and 
Fig. 6, B 2 and Bo are independent of 7(= ^3), and only B\ depends on j(= 3 ). Three 
different lines for B\ correspond to the relevant values of 7(= ^3), respectively. It is clearly 
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7=65 - 
7 = 85 " 




23456789 10 

Nc 

FIG. 5. The branching ratio (in units of 10~ 4 ) versus the effective number of color, N c , for 
B° — > ir~ X decay. B tot (= B) has been calculated using \V u b\ = 0.004 and is denoted by the bold 
solid line. The solid, dotted, and dashed lines correspond to 7 = 65°, 85°, 110°, respectively. For 
light quark masses, m u = 5 MeV and = 7 MeV have been used. 



shown that £> 2 is dominant. A representative value of B for \V u b\ = 0.004 and 7(= 3 ) = 85° 
is shown as the bold solid line in the figure. The value of B does not vary much as iV c varies. 

Similarly, Figure 6 shows the BR of 5° — > ix~X as a function of 7(= 3 ) for three 
different values of N c = 2, 5, 10. The solid line corresponds to the case N c = 2, and the 
dotted line and the dashed line are for N c = 5 and N c = 10, respectively. The £> 2 , B\ and 
Bo increase as iV c increases. However, since the dominant term B2 does not change much 
when N c is varied, the BR does not change much either. In Fig. 6, a representative value of 
B for \ V u b\ = 0.004 and N c = 5 is shown as the bold solid line. 

Finally we summarize our result in Fig. 7. The BR of B° — > n~X is presented as a 
function of \V u b\. For light quark masses, we use m u = (1.5 — 5) MeV and m d = (3 — 9) 
MeV. We also vary the value of iV c and 7(= 3 ) in a reasonable range: from N c = 2 to 10, 
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TABLE II. The branching ratio of B -> tt _ X for fixed N c = 3. Here B denotes the BR 
calculated for \V ub \ = 0.004. 



7 


B x 10 4 


B 2 x 10 4 Bi x 10 4 


B x 10 4 


60 


1.18 


0.98 


0.18 


0.019 


70 


1.16 


0.98 


0.16 


0.019 


80 


1.14 


0.98 


0.14 


0.019 


90 


1.11 


0.98 


0.11 


0.019 


100 


1.07 


0.98 


0.072 


0.019 


110 


1.04 


0.98 


0.037 


0.019 


TABLE III. The branching 


ratio of B° — > ir~X for fixed 7 = 65° and iV c = 3. 


Here B denotes 


the BR calculated for \V ub \ = 0.004 and 2.32GeV < E n < 2.56GeV. 




7 = 


65° 




N c = 3 




N c B x 10 


4 


7 


B x 10 4 


2 0.95 




60 


1.06 


3 1.05 




80 


1.04 


4 1.11 




100 


1.01 


5 1.14 




110 


0.99 



and from 7(= 3 ) = 60° to 110°. The solid and the dotted lines correspond to the smallest 
and the largest value of B in the given parameter space, respectively. For the given \V ub \, 
the BR is estimated with a relatively small error (< 20%), as can be seen. Reversely, for 
the given (i.e., experimentally measured) BR, the value of \V ub \ can be determined with a 
reasonably small error (~ 15%). (Of course, since in practice the BR would be measured 
with some errors, \V u b\ could be determined with larger error: e.g., for B = (1.0±0.1) x 1CT 4 , 
our result suggests \V ub \ =~ (3.7 ± 0.59) x 10~ 3 .) 

In order to use the decay process 5° — > n~X, one may need to consider the B° — 5° 
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Nc-5 " 
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FIG. 6. The branching ratio (in units of 10 -4 ) versus the CP phase, 7, for B — ► ir~X decay. 
Btot{= B) has been calculated using \V u b\ = 0.004 and is denoted by the bold solid line. The solid, 
dotted, and dashed lines correspond to iV c = 2, 5, 10, respectively. m u = 5 MeV and = 7 MeV 
have been used. 



mixing effect: 5° — > B° — > n~X. The neutral B° has about 18% probability of decaying as 
the opposite flavor B° [M]. Thus, including the B° — B° mixing effect, the decay rate T 
for 5° decay to n~X can be expressed as 

T = (0.82) • T(B° -> n~X) + (0.18) • T(B° -> tt - X) , (23) 

where r(5°(5°) — > 7r _ X) denotes the decay rate for B~°(B°) decay directly to tt~X. The 
BR of B° decay directly to n~X is about 90% of the BR of B° decay directly to tt~X with 
the energy cutQ, 2.32GeV < E w < 2.56GeV, as can be seen from Table III. Even though the 



3 As mentioned in Sec. II, the charged pion in the decay mode B° — > ir~ X contains the spectator 
quark d, and this process is basically a three-body decay B° — ► ir~ud (s). Therefore, in order to 
remove this large s quark contribution, one needs to make such a large energy cut (see Fig. 4). 
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2.0 
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2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 4.6 4,8 

|Vub| (xlO 3 ) 

FIG. 7. The branching ratio (in units of 10 -4 ) versus | T^t. I for B° — > 7r~X decay. The solid 
and the dotted lines correspond to the smallest and the largest value of B in the given parameter 
space, respectively. For light quark masses, m u = (1.5 — 5) MeV and = (3 — 9) MeV have been 
used. 

theoretical estimate for the mode B° — > tt~X would include a somewhat larger uncertainty, 
which mainly arises from the relevant hadronic form factor, the total error of T would not 
increase much. For example, if the estimate of T(B° — > tt~X) in Eq. (|2BJ) includes an error 
of 30%, then its actual contribution to the final error of Tq is less than 5%. Therefore, 
even after considering the effect from the B° — 5° mixing, our result holds with reasonable 
accuracy. 

IV. CONCLUSIONS 

We have studied semi-inclusive charmless decays of B— mesons to nX in the final state, 
such as B° — > n^^X, B° — > 7r ± ^X, 5 ± — > tx ± ^X 1 where X does not contain a charm 
(anti)quark. Among these B — > ttX decays, we have found that the mode 5° — > tt~X 
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(B° — > n + X) is particularly interesting and can be used to determine of the CKM matrix 
element \V u b\ m phenomenological studies. 

In B° — > 7T~X decay, the charged pion in the final state can be produced via W— boson 
emission at tree level and is expected to be energetic (E w ~ m B /2). Thus, the energetic 
charged pion in the final state can be a characteristic signal for this mode. This process 
is basically a two-body decay process of b — > ix~u. As a result, in this mode, the model- 
dependence does not appear to be severe. 

We have calculated the BR of 5° — > 7r~X and presented it as a function of \V u b\. It is 
expected that its BR is an order of 10~ 4 . (In this analysis, higher-order QCD corrections 
have not been considered; instead we analyzed within the QCD improved general factoriza- 
tion framework. So, a further study on this process would be very interesting.) We have 
also estimated the possible uncertainty due to B° — 5° mixing effects via the decay chain 
B° — > B° — > n~X. Other theoretical uncertainties, such as those arising from the WC's and 
the CKM elements, could affect our results in some extent. However, as soon as the rele- 
vant results from the experiments become available, one can use them to reduce theoretical 

uncertainties in turn. Thus, in the viewpoint of phenomenological studies, our results can 

— o 

be used to determine \V u b\ with reasonable accuracy by measuring the BR of B — > tx~X. 
Therefore, the process B U -> TT-X (B° -> 7T + X) can play an important role in measuring 
| V ub | at B— factories. 
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